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Abstract 

We have investigated the possibility of inferring peculiar velocities 
for clusters of galaxies from the Doppler shift of scattered cosmic mi¬ 
crowave background (CMB) photons. We hud that if the core radius 
of the gas distribution or the beam size of the instrument is larger 
than 3-7 arcminutes, then the maximum attainable signal-to-noise ra¬ 
tio is determined by confusion with primary fluctuations. For smaller 
angular scales, “cosmic confusion” is less important and instrumental 
noise and/or foreground emission will be the limiting factor. For a 
cluster with the optical depth of the Coma cluster and for an optimal 
Ritering technique, typical one-sigma errors span the wide range from 
400 to 1600 kms“^, depending on the cosmological model, the resolu¬ 
tion of the instrument and the core radius of the cluster. The results 
have important implications for the design of future high-resolution 
surveys of the CMB. Individual peculiar velocities will be measurable 
only for a few fast moving clusters at intermediate redshift unless cos¬ 
mic fluctuations are smaller than most standard cosmological scenarios 
predict. However, a reliable measurement of bulk velocities of en¬ 
sembles of X-ray bright clusters will be possible on very large scales 
(100 - 500/j-i Mpc). 



1 Introduction 

In a seminal paper, Sunyaev & Zeldovich (1970) discussed the important ef¬ 
fect of Compton scattering of cosmic microwave background (CMB) photons 
by hot electrons. The ionized intra-cluster gas of galaxy clusters leads to a 
Compton distortion of the CMB behind the cluster on scales of a few arcmin- 
utes, and this effect is now generally referred to as the SZ effect (Sunyaev & 
Zeldovich 1972, 1980a). After a twenty year long history of only tentative 
detections, the number of convincing observations is now rapidly increas¬ 
ing (Birkinshaw, Hughes & Arnaud 1991; Klein et al. 1991; Grainge et al. 
1993; Jones et al. 1993; Birkinshaw & Hughes 1994; Wilbanks et al. 1994). 
Detector technology and our understanding of the numerous noise sources 
is developing fast, and it should soon become possible to observe some of 
the “secondary” effects of this scattering process which were predicted by 
Sunyaev & Zeldovich (1972, 1980b). The most prominent of these is an ad¬ 
ditional Doppler shift of the CMB signal due to the radial peculiar motion 
of the cluster with respect to the CMB rest frame (see Rephaeli & Lahav 
(1991) for a recent comprehensive discussion). The amplitude of the corre¬ 
sponding CMB distortion is expected to be a factor 10-20 smaller than that 
for the “ordinary” SZ effect. Following Rephaeli & Lahav (1991), we will re¬ 
fer to the two effects as the thermal and the kinematic SZ effect, respectively. 
No convincing detection of the kinematic effect has yet been reported. The 
best hope for a detection comes from proposed second generation satellite 
missions, which would reduce pixel noise levels by more than a factor of ten 
compared to COBE and provide much higher angular resolutions, perhaps 
ten arcminutes or better — see e.g. Mandolesi et al. (1995). In addition, a 
number of ground-based experiments with suitable angular resolution and 
sufficient sensitivity are currently being built or planned for the near future. 

The potential pay-off from a reliable measurement of cluster peculiar 
velocities with respect to the CMB frame would be enormous. So far, 
not much is known about bulk velocities on very large scales (redshifts 
> 10000 km s“^). The measurement of the bulk velocity of a volume-limited 
sample of 119 Abell clusters out to a distance of 15000 km s“^ gave a rather 
high value of 700 — 800kms“^ (Lauer & Postman 1994; Postman & Lauer 
1995; Colless 1995, see Strauss & Willick 1995 for a recent general review 
on peculiar velocities). This is considerably larger than the approximately 
400 —500kms“^ expected in most cosmological scenarios currently discussed 
(Feldman & Watkins 1994; Strauss et al. 1995). An independent measure¬ 
ment method on even larger scales should clarify the situation. 
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As we will discuss below, on top of technical difficulties and the com¬ 
plications caused by galactic and atmospheric foreground emission, there is 
an additional important impediment to a measurement of the kinematic SZ 
effect which so far has received little or no attention: the possible confusion 
of a Doppler shift of the CMB due to the peculiar velocity of the cluster 
with primary fluctuations of the CMB on the same scales. This will prob¬ 
ably turn out to be a fundamental limitation for a determination of cluster 
peculiar velocities by the kinematic SZ effect with beam sizes realizable with 
space missions. In this paper we will review the basic principles of the effect, 
discuss briefly the numerous problems involved with detecting it, give a de¬ 
tailed assessment of the possible confusion with primary CMB fluctuations 
and discuss what results future measurements may be able to produce. 


2 Measuring the kinematic SZ effect 


2.1 Its strength and spectral signature 


Peculiar velocities of the hot intra-cluster gas lead to a Doppler shift of 
the scattered photons which is proportional to the product of the radial 
peculiar velocity Upec and the electron density integrated along the line 
of sight through the cluster. For small optical depths, the distortion from 
a blackbody spectrum is small and the relative change in intensity of the 
CMB is given by 
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where T is the temperature of the microwave background and ctt is the 
Thompson cross section (Sunyaev & Zeldovich 1972a,b). In the Rayleigh- 
Jeans limit a; ^ 1, this is approximately 
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Equation (1) is very similar to the one for the thermal effect, where Upec/c is 
replaced by kT^lm^c^ and the spectral dependence of the relative intensity 
change has an extra factor [x coth(a;/2) — 4]. 

For a typical cluster, the kinematic SZ effect at the cluster center is of 
order 
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where is the electron density in the core, is the core radius and we 
have scaled to the values reported by Briel et ah 1991 for the Coma cluster 
(assuming a distance of 140 Mpc). 

2.2 Noise sources 

The actual observed signal will, however, be a convolution with the beam 
pattern -B(x) of the experiment, and will contain a considerable noise com¬ 
ponent. 


f^(x)j = y 5(x - x')(5(x') iV,ky(x'))d^x' iVdet(x), (4) 

where the noise-free unconvolved signal ^(x) will be given by equation (1). 
iVgky and iVdet denote the noise in the sky and in the detector, respectively. 
The area of integration varies considerably between different experiments 
and can be rather complicated. Most ground-based observations use drift 
scans and beam switching, while the proposed space-borne missions will 
scan the whole sky in a complicated pattern. There is a wide variety of 
noise sources which make a measurement of the kinematic SZ effect rather 
difficult. The most important are the following: 

• Instrumental noise 

• Diffuse foreground emission (atmospheric, galactic, extragalactic) 

• Radio sources 

• The thermal SZ effect of the cluster 

• Bulk motion of the gas within the cluster 

• Primary CMB fluctuations 

In the remainder of this section, we will give a very brief review of the 
hrst hve of the above noise sources (see Fisher & Lange 1993 for a more 
comprehensive discussion). We will return to the last item, the possible 
confusion with primary fluctuations, in Section 3. 

2.2.1 Instrumental noise 

The instrumental noise is mainly due to the electronics of the receivers. It is 
usually assumed to be independent from pixel to pixel. This corresponds to 
a spatial power spectrum equivalent to that of white noise — see Tegmark 
& Efstathiou (1995), hereafter “TE95”, for a detailed discussion. For a 
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space-borne mission, the noise level currently achievable is of the order of 
5 (0/10')“^//K at 200 GHz (Mandolesi 1995). 

2.2.2 Diffuse foreground emission 

Atmospheric emission due to water vapour in the atmosphere will be at a 
level of about 50 — 100 //K (Fischer & Lange 1993) at a frequency of 250 
GHz and below. At higher frequencies, it increases rapidly. This will make 
measurements of the kinematic SZ effect from the ground very difficult. 

Outside the Earth’s atmosphere, galactic synchrotron radiation and brems- 
strahlung are the main sources of diffuse foreground emission at low fre¬ 
quencies (below the peak in the CMB spectrum), while at high frequencies, 
infrared emission from galactic dust and (young) galaxies is dominant. As 
we will discuss further on, the most favourable frequency for observing the 
kinematic SZ effect is somewhat above the peak in the CMB spectrum, at 
approximately 220 GHz. Here only the latter two effects potentially pose 
problems. 

Our best information about the small-scale properties of the galactic 
dust comes from the 100, 60 and 25 // IRAS maps. An extrapolation from 
these very high frequencies involves considerable uncertainties. However, 
the measurements of the kinematic SZ effect will only be sensitive to varia¬ 
tion of the galactic dust emission on the rather small angular scales of ten 
arcminutes and below. Recent analyses of the IRAS maps (Low & Cutri 
1994, TE95) indicate a spatial power spectrum for the galactic dust emis¬ 
sion for which the power per decade in scale falls off rapidly towards smaller 
scales. A rough estimate of the expected dust contamination on 10’ scales 
around 200 GHz is 10 //K in fairly clean parts of the sky (TE95). Since this 
contamination can be hltered out to some extent by using multi-frequency 
data, it is thus likely that dust will cause slightly less of a problem than 
instrumental noise. 

The far infrared background emission due to the combined emission from 
(young) galaxies is hard to predict, as little is known about the properties 
of galaxies in the relevant wavelength range. Furthermore the evolution of 
number density, luminosity and spectrum of these objects towards higher 
redshift is highly uncertain. Theoretical estimates vary between a few and a 
few tens of ^uK, comparable to or slightly above the expected instrumental 
noise (Guiderdoni et al. 1995). 
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2.2.3 Radio sources 


Discrete radio sources are usually one of the major difficulties for measure¬ 
ments of the Sunyaev-Zeldovich effect from the ground. However, these 
measurements are generally carried out at lower frequencies. The spectrum 
of the CMB rises rapidly (///^ ^ below of the peak of the spectrum at 160 
GHz), much faster than typical radio source spectra. Radio source contam¬ 
ination will therefore be much less of a problem at the optimal frequency of 
218 GHz. 

2.2.4 The thermal SZ effect 

The thermal SZ effect is generally considerably larger than the kinematic SZ 
effect. The exact factor of course depends on the density-weighted electron 
temperature and the peculiar velocity of the cluster. In the Raleigh-Jeans 
limit, the corresponding ratio of the observed change in brightness temper¬ 
ature is approximately a factor 20 (re/10®ld)(npec/500 kms“^)“^. However, 
as discussed in Section 2.1, the two effects are distinguishable by their differ¬ 
ent spectra. In the case of the thermal SZ effect, the change in the observed 
brightness temperature is a Compton distortion with a decrease at low fre¬ 
quencies, an increase at high frequencies and a zero-crossing at ~ 218 GHz. 
The kinematic SZ effect, in contrast, has the same sign for all frequencies. 

2.2.5 Bulk motions within the cluster 

The bulk motion within the cluster (including a possible rotation) of the 
hot gas is expected to be small if the cluster is fully relaxed and the gas 
has attained hydrostatic equilibrium. However, evidence is increasing that 
galaxy clusters spend only a fraction of their lifetime in such an equilibrium 
state due to frequent merger events (Navarro, Frenk & White 1995). Internal 
bulk motions of the gas could therefore be frequent and of considerable 
amplitude. While these motions and the corresponding Doppler effect will 
be less relevant for poorly resolved or unresolved objects, they might be an 
important source of noise when hltering techniques are applied for extended 
objects. We will give a short discussion of this effect in Section 4.3. 

2.2.6 Primary fluctuations 

Observations of the DMR experiment on board the CODE satellite have 
established the existence of temperature fluctuations of the CMB on scales 
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of a few degrees and above and determined their amplitude to be of order 
(AT/T) ~ 10“® (Smoot et al. 1992). In the last two years, this result has 
been conhrmed by a number of ground- and balloon-based experiments sen¬ 
sitive to temperature fluctuations on angular scales from a few degrees down 
to a few arcminutes. The scatter between the different measurements is still 
about a factor of two, but the measured amplitudes are consistent with 
the assumption that the fluctuations are generated by primordial potential 
fluctuations extending over all scales. In standard cosmological scenarios as¬ 
suming a flat (Harrison-Zeldovich) power spectrum of potential fluctuations, 
the amplitude of temperature fluctuations should rise towards smaller scales 
and reach a maximum around 10 arcminutes (a scale set by the horizon size 
at the time when radiation and matter-density were equal). At even smaller 
scales a dramatic drop should occur due to diffusive damping (White, Scott 
& Silk 1994 for a recent general review on the CMB and Hu, Sugiyama & 
Silk 1995 for a description of the physics responsible for the Doppler peaks). 
Unfortunately, the spectral distribution of these intensity fluctuations is ex¬ 
actly the same as that of the kinematic SZ effect, be., they correspond to a 
blackbody spectrum whose temperature varies slightly with position in the 
sky. For scales larger than the diffusive damping scale, cosmic fluctuations 
will therefore be the dominant noise source. 

2.3 Isolating the kinematic SZ effect 

The presence of the numerous noise sources described above make a multi¬ 
wavelength measurement absolutely essential to get a handle on possible 
systematic errors. In principle, modeling of noise components with different 
spectral signature should allow a discrimination of as many components as 
there are observed frequencies. In practice, when faced with both detec¬ 
tor noise and considerable uncertainties in such modeling, a more robust 
approach is to choose a suitable frequency and estimate the impact of the 
dominant sources of contamination at that frequency. The simplest strategy 
is probably the use of a channel centered around the cross-over frequency 
of the thermal SZ effect (218 GHz), so that negative and positive contri¬ 
butions from lower and higher frequencies effectively cancel (Rephaeli & 
Lahav 1991). This channel would still be in the favourable frequency range 
where the diffuse foreground contamination is small. As an additional ben- 
eht, the radio source contamination will be much smaller than at the lower 
frequencies generally used so far for SZ measurements from the ground (cf. 
Willbanks et al. 1994). 
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3 Optimal filtering to minimize confusion with “cos¬ 
mic” fluctuations 

3.1 Estimating the noise due to cosmic fluctuations 

3.1.1 The estimate for the peculiar velocity and its error 

To get an estimate Eest for the radial peculiar velocity of a cluster, we need 
some information about the gas distribution in the cluster and the corre¬ 
sponding optical depth (which today would be obtained from X-ray observa¬ 
tions, but in the future might also come from observation of the thermal SZ 
effect). If the cluster is resolved, we can choose an appropriately normalized 
hlter function and choose our estimate to be 



The hlter function should be normalized so that 

Upec = / S(x^) -(/’(x) d?x (6) 

Jf 

is the true density weighted peculiar velocity of the cluster — '^(x) again 
denotes the noise-free signal. As shown in the appendix, the mean squared 
error of the estimator is given by 

((Kst - Vc)^> = ^ X 

if the CMB “noise” properties are translationally invariant and can be de¬ 
scribed by a power spectrum T’(k). Here and throughout, hats denote 
Fourier transforms. 


3.1.2 The assumed cluster profile and a first estimate of the ex¬ 
pected noise level 


To demonstrate the most important effects, we will use axisymmetric /3- 
models to model the distribution of the cluster gas; 


rae(r) = UeO 
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where rieo is the central electron density and is the core radius (Cavaliere 
& Fusco-Femiano 1976, 1978). With 0.5 < /3 < 1, these models are generally 
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a good fit to the electron density inferred from X-ray observations (Jones & 
Forman 1984). The optical depth of the cluster is then given by 


r(x) = Tig Gt dl = UeO (Tt J’c 


r(3/3/2- 1/2) 
r(3/3/2) 


1 + 


• (9) 


Typical central optical depths for luminous X-ray clusters are of order 0.002 — 
0.02. We now calculate the expected signal-to-noise as a function of the core 
radius of the cluster using a “naive” filter function ^(x) oc t'(x) appropri¬ 
ately normalized to get a first estimate of the expected noise level. t'(x) 
denotes the convolution of the cluster profile with the beam profile. 

In Figure 1, the corresponding la error is shown for three different cos¬ 
mological models. The following points should be noted here: 

• The overall noise level is of order 250 //K, corresponding to about 
1500 (r/0.02)“^ km s“^ (rather independent of the core radius for large 
core radii). This is a factor 10-50 higher than the other noise sources 
discussed in Section 2, and even for this optimistic optical depth, the 
expected la error in the peculiar velocity measurement is rather large 
compared to typical cluster velocities. 

• Beam dilution sets in when the core radius of the cluster becomes 
smaller than the resolution of the instrument. The signal-to-noise then 
scales as 0“^. 

• The overall noise level depends only very weakly on the cosmological 
model. 


3.2 Optimal filtering 

3.2.1 The general shape of the filter function 


A glance at equation (7) shows that the expected signal-to-noise depends 
on the filter function chosen. If we know the noise power spectrum well, we 
can optimize the filter function by minimizing the variance of our estimator. 
This is a simple variational problem, whose solution is given in Appendix A. 
We find that the Fourier transform of the optimal filter function 'ip is given 

by 
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where r' is the Fourier transform of r'. The corresponding signal-to-noise is 
given by 



For Gaussian noise, the interpretation is simple, as the error in the estimate 
will then also be Gaussian distributed. Figure 2 shows the prohle of typical 
hlter functions. It has alternating positive and negative ’’sidelobes” which 
effectively cancel the long-wavelength cosmic fluctuations. The cluster pro¬ 
hle is shown for comparison (both in arbitrary units). 

3.2.2 Axisymmetric filter 

To demonstrate the main effects of optimal hltering, we again use axisym¬ 
metric /3-prohles to characterize the radial prohle of the gas distribution. 
Figures 3a-c show the effect of optimal hltering on the expected la error of 
the estimate for the peculiar velocity as a function of the core radius of the 
cluster. The overall noise level is reduced by approximately a factor of two 
on all scales. This is due to the broad-band character of the cosmic huctu- 
ations, with a substantial contribution from huctuations with wavelengths 
much larger than the characteristic scale of the cluster — wavelengths which 
can efficiently be hltered out. However, the optimal hltering technique be¬ 
comes even more efficient if the diffusive damping scale, below which the 
cosmic huctuations are essentially zero, is above the characteristic scale of 
the cluster. The noise level can then be reduced to the pixel noise of the 
instrument if the instrumental resolution is sufficient. Figures 3a-c show the 
inhuence of the instrumental resolution, the pixel noise and the details of 
the cosmological model on the signal-to-noise ratio. The achieved signal-to- 
noise depends very sensitively on the size of the core radius of the cluster 
relative to the resolution of the instrument and the diffusive damping scale. 

3.2.3 Non-axisymmetric filter 

If a high quality X-ray map of the cluster is available, we can do better 
than htting a simple /3-proffie to the surface brightness distribution and 
constructing an axisymmetric hlter function. To demonstrate this, we have 
used a surface brightness map of Coma (kindly provided by Hans Bohringer) 
and used the full two-dimensional spatial information to construct a non- 
axisymmetric spatial hlter function. In Figure 4, the results are compared 
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to axisymmetric filtering for different beam sizes and different values of the 
pixel noise. The cosmic noise is efficiently filtered out even for core radii 
comparable to or larger than the scale of the Doppler peaks. The achievable 
signal-to-noise now depends on the level of instrumental noise for all scales 
(as well as on the resolution of the instrument). However, as discussed in the 
next section, a spatial filter function which is optimized using such detailed 
information about the cluster gas distribution is more sensitive to errors in 
the assumed distribution. 

3.2.4 Possible complications 

To use a non-axisymmetric optimal filter function, we need good knowledge 
of the spatial distribution of the line-of-sight integral of the radial velocity 
times the density f t>pec(x) ra(x)d/. From the surface brightness distribution 
of a high-quality X-ray map we will get an estimate of the optical depth 
distribution. However, when constructing a spatial filter function, we will 
have to neglect the errors due to the unknown internal bulk motions, a pos¬ 
sible clumping of the gas, the noise in the X-ray map and so on. To get at 
least a rough feeling for the effect of such errors in the assumed distribution 
of the signal, we have used clusters which were simulated with a Smoothed 
Particle Hydrodynamics code (kindly provided by Matthias Steinmetz) and 
were assumed to move with a peculiar velocity of 1000 km s“^. The simula¬ 
tions are similar to those described in Navarro, Frenk & White (1995) and 
show internal bulk motions of order 500kms“^. We found that the esti¬ 
mated peculiar velocities were systematically smaller than the true optical 
depth-weighted mean values. The mean value of fest/Drue was 0.89 with a 
standard deviation of 0.18 for a sample of 12 cluster seen from three different 
directions. This should give an indication of the size of the effect of errors 
in the assumed signal distribution if non-axisymmetric filtering is applied. 
For axisymmetric filter functions, the errors were considerably smaller. 

4 Determining cluster peculiar velocities 

4.1 Using X-ray data to estimate the expected errors 

The key parameters for a measurement of the kinematic SZ effect are the 
optical depth in the core, the core radius and the peculiar velocity of the 
cluster. In Appendix B, the luminosity-temperature relationship of galaxy 
clusters is used to derive the typical scaling of optical depth and core radius 
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with X-ray luminosity, r ~ -^x ~ These scaling laws can be 

used to estimate optical depth and core radius relative to those of a reference 
cluster from the ratio of the X-ray luminosities and redshifts. In Figure 5, 
such estimates are shown for the XBACs cluster sample (X-ray brightest 
Abell-type clusters exceeding a flux limit of 5.5 X 10“^^ erg cm“^ s“^ in the 
ROSAT band; Ebeling 1993, Ebeling et ah 1995). With the exception of a 
few well-studied clusters, the core radius and the central optical depth of a 
cluster are poorly determined quantities. We have therefore chosen Coma 
as a reference cluster for which we assume an optical depth r = 0.01 and 
a core radius of 10 arcminutes (Briel et al. 1992; see Herbig, Readhead & 
Lawrence 1992 for a measurement of the thermal SZ effect). However, one 
should note that these estimates will only be correct in a statistical sense, as 
there is a rather large spread in the luminosity-temperature relation of X-ray 
clusters reflecting the fact that there is no universal gas mass distribution for 
a cluster of given luminosity and redshift. Especially, cooling flow clusters 
will deviate from these scaling laws and will generally have smaller and 
denser cores (Edge & Stewart 1991; Fabian et al. 1994). Using Coma as a 
reference cluster should therefore yield conservative error estimates. 

As expected, the angular extent of the core radii are generally smaller 
than for Coma, mainly due to the larger distances. Typical values lie be¬ 
tween 1 and 3 arcminutes, corresponding to a typical redshift around x = 0.1. 
The optical depth lies between one hfth and twice that of Coma. We can 
now use the derived scaling laws to compute the expected la error of a 
measurement of the kinematic SZ effect for the same cluster sample as 
above. In Figure 6, these are shown for our standard cosmological model 
(fljQj = l,flj = 0.06) and two different beam sizes (4 and 8 arcminutes). 
Here we have assumed that an X-ray map will exist of sufficient quality to 
ht a /3-model to the X-ray surface brightness, so that we can apply the hl- 
tering technique with an axisymmetric hlter. As discussed above, there is 
a characteristic angular scale, above which confusion with primary fluctu¬ 
ations becomes important and which is about 5 arcminutes for the special 
cosmological model chosen here. The la errors drop considerably for clus¬ 
ters with smaller core radii if the beam size is sufficiently small to “avoid 
the Doppler peaks”. Even though this is difficult, it might be possible for a 
future space-borne mission. For a diffraction limited instrument with aper¬ 
ture B observing at frequency //, the maximum resolution is approximately 
5'(D/1 m)-i(///200 GHz)-U 


11 



4.2 Individual clusters 

For standard cosmological scenarios, it seems impossible to obtain more than 
upper limits of order 1000 km s“^ for the peculiar velocity of clusters with 
core radii larger than 3-7 arcminutes, even if the gas distribution is known 
well enough to apply non-axisymmetric hltering. This probably precludes 
a useful determination of the peculiar velocity of a nearby extended X-ray 
luminous cluster like Coma. However, chances are better for nearby cooling- 
flow clusters which have smaller core radii and a larger central optical depth 
(Edge & Stewart 1991). Generally, prospects are best for the most X-ray 
luminous clusters at somewhat higher redshifts, with core radii around 3 
arcminutes. Figure 5 shows that there are approximately 30 such suitable 
clusters with a core radius around 3 arcminutes and r 0.01. The mini¬ 
mum 1(7 error should here be around ~ 300 — 700 km s“^. This might just be 
sufficient if the peculiar velocities are as high as indicated by recent peculiar 
velocity studies of Abell clusters out to a redshift of 15000 km s“^ (Lauer 
& Postman 1994; Postman & Lauer 1995; Colless 1995). For a diffraction- 
limited instrument with an aperture of about 1 m, the measurement of the 
kinematic SZ effect for clusters with smaller core radii will suffer consider¬ 
ably from beam dilution effects. The minimum measurable peculiar velocity 
therefore increases again towards very small core radii. 

Figure 7 shows the distribution of la errors of the cluster sample for 
different cosmological models and a beam size of 4 arcminutes. The errors 
generally increase if fftot is lowered and decrease if ffbar is lowered. 

4.3 Statistical measurement of bulk velocities on large scales 

An obvious way to improve the signal-to-noise ratio is averaging the signal 
of ensembles of clusters. To get a feeling for what could be achieved, we 
have calculated the expected error for a maximum likelihood solution of 
the bulk flow velocity with respect to the microwave background, again for 
the XBACs cluster sample (Ebeling 1993). Neglecting small-scale peculiar 
velocities, the maximum likelihood solution for the bulk flow velocity U of a 
sample of a clusters at positions r^ with measured radial peculiar velocities 
Vg is given by 


r-'. = 4y S' 
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where r = r/|r| and the are the rms peculiar velocity errors (Kaiser 1988). 
The covariance matrix of the maximum likelihood estimates is 

([/,[/,) = A-d. (13) 

Table 1 shows the rms errors for the x, y and z components of U in galactic 
coordinates for a set of cosmological models and two beam sizes. The number 
of clusters in the sample is 244 and the mean redshift of the sample is 
2 ; = 0.088. The values indicate that a reliable measurement of bulk velocities 
should be possible at redshifts, cz ~ 10000 — 50000 km s“^ corresponding to 
scales of 100 — 500h“^ Mpc. 

5 Conclusions 

In view of the severe obstacle of atmospheric noise facing ground-based in¬ 
struments, a space-borne mission seems to be the most promising method 
to measure the kinematic SZ effect. As pointed out previously, mapping the 
sky at a frequency of 218 GHZ, where the thermal SZ effect vanishes, should 
be the optimal strategy. The dominant noise source will here be confusion 
with primary fluctuations. The expected noise level is of order or larger 
than the expected signal, and it is therefore essential to use the knowledge 
of the CMB “noise” properties and the gas distributions of the individual 
clusters (which can be obtained by the mission itself and from X-ray ob¬ 
servations, respectively). This knowledge makes it possible to analyze the 
CMB maps with a spatial hlter optimized for individual clusters. An im¬ 
provement in signal-to-noise by a factor of two is easily achievable, and even 
a factor of 10 is possible if the gas mass distribution is well known from a 
high-quality X-ray map. The hnal signal-to-noise ratio depends crucially on 
the cosmological model and the angular resolution of the instrument, and it 
is not currently clear whether a meaningful peculiar velocity measurement 
for individual clusters will be possible. Prime candidates are X-ray lumi¬ 
nous clusters at intermediate redshift with core radii just below the Doppler 
peak scale. For a favourable but still rather standard cosmological scenario 
(standard CDM with low baryon fraction) and a good angular resolution (4’ 
FWHM), the peculiar velocity of as many as 30 individual clusters might be 
determined accurately. Even if this is impractical, it should still be possible 
to determine the bulk motion of an ensemble of 200 X-ray luminous clusters 
at redshifts ^ 10000kms“^ with an accuracy of order 200kms“^. 
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Appendix A: Optimal filtering 

In this Appendix, we derive the optimal hltering scheme which was used in 
the main part of the paper. Suppose we want to extract the amplitude A 
of a signal with a known spatial distribution t(x) from a measured signal 
Il(x) which is contaminated by noise iV(x). Writing 

I)(x) = A r(x) + iV(x), (14) 

the application in this paper corresponds to t(x) being the known optical 
depth prohle of some cluster, A being proportional to the radial velocity fpec, 
and the random held N being the nuisance contribution from the CMB and 
pixel noise. Let us make no other assumptions about N than that it has 
zero mean and that its statistical properties are independent of position. 
Letting hats denote Fourier transforms, this implies that 

(iV(x)> = 0, (15) 

(iV(k)*iV(k)) = (27r)2^(k'- k)P(k), (16) 

for some power spectrum P(k). Note that we are not making any assump¬ 
tions about N being Gaussian. The most general estimator of A that is 
linear in the data can clearly be written as 

Aest = [ il){x)D{x)£x (17) 


for some weight function ijj. Using equations (15) and (16), we readily obtain 


b 


(Aest - A) = j V’(x)r(x)d^a; - 1, 


((Ae.t - Af) = h^A 


(2vr) 


V’(k) P(k)£k, 


(18) 

(19) 


were b, the bias, gives the average error and a, the noise, gives the r.m.s. 
deviation of the estimate from the true value. The goal of this Appendix is 
to hud the best estimator Aest, be., the optimal choice of ijj. Now what do 
we mean by an estimator being good! We will use the following criteria: 

• Require it to be unbiased {b = 0) 

• Minimize the noise a 
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The requirement of no bias clearly just fixes the normalization of whereas 
its shape will be determined by the second condition, the requirement that 
the average error be minimized. The problem of minimizing subject 
to the constraint that 5 = 0 is readily solved by introducing a Lagrange 
multiplier A and hnding the function 'ip that minimizes T = + Xb. Using 

Parseval’s theorem to re-express b in terms of we thus seek the function 
p} that minimizes 

L = ^ I my [t^(k)P(k) + Af(k)] £k - A, (20) 


Taking the functional derivative with respect to y and setting the result 
equal to zero, we thus obtain the simple result ^(k) oc f(k)/P(k). This 
has a simple intuitive interpretation: the optimal hlter y is such that it 
is maximally sensitive to those wave numbers k where f is large and P is 
small. In other words, it focuses on those Fourier modes where the cluster 
prohle stands out the most above the noisy background. Normalizing this 
so that b = 0, we get 


^(k) 


1 f ^ 

.myj p(k') J p(k)’ 


( 21 ) 


our desired result. Substituting this back into equation (19), we hnd that 
the best signal-to-noise ratio A/a that can be attained is 


A 

a 



l^(k)P 

P(k) 


cPk 



( 22 ) 


To apply these formulas in real-world situations, the effects of pixel noise and 
beam smoothing are simply absorbed into P{k) as described in T95. The 
quantities of equation (14) all become convolved with the beam function P, 
and N receives an additional contribution from pixel noise, leaving 


P'(x) = zl r'(x)-|-iV'(x). (23) 

Here r' denotes r convolved with the beam as before (be., f' = Pr), D' 
denotes the noisy beam-smoothed signal that we actually observe, and as 
shown in T95, N' has the power spectrum 

P'(k) = P,,,b(k)|P(k)p + ^. (24) 
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Here ra/47r is the number of pixels per steradian, and a„ denotes the r.m.s. 
noise level per pixel. Substituting this into the solution we found above, we 
thus obtain the optimal hlter function ^(k) oc f'(k)/P'(k), i.e., 

_fOO_ 

Thus the best signal-to-noise ratio A/a that can be attained is 


^(k) = 


|f(k')|^ 


(2vr) 


P 


4.7ra‘z 


-d^k' 


cmb v 


n\B(k'M^ 


A 

a 


1 

(27r)2 


l^(k)P 


1 1/2 


-Fcmb(k) + 


-d'^k 


n|B(k)P 


A. 


If the instrument has a Gaussian beam, then 


(26) 


P(k) exp [—6^ (/ + 1) /] , (27) 

where 9 = FWHM/2-\/ln 2 0.6FWHM and FWHM denotes the full-width- 
half-max. 


Appendix B: Using scaling relations to estimate the 
optical depth and core radius of a cluster 

It is generally accepted that the properties of clusters of galaxies obey certain 
scaling relations, even though the scatter is large (Edge & Stewart 1991). 
The scaling of the X-ray luminosity with the temperature of the emitting 
gas is the best established of these correlations. Lx ~ T“ with a ~ 3. 

For bremsstrahlung emission from a cluster of temperature T with core 
radius and core density ric, the X-ray luminosity should scale as 

Lx ~ nl (28) 

The temperature of the gas should be determined by the depth of the poten¬ 
tial well of the gravitationally dominant dark matter halo in which the gas 
is believed to reside. If the gas traces the dark matter (as indicated by nu¬ 
merical simulations), the temperature of the gas should scale as ~ r^, 

which translates into 

Lx^T-^l'^r;\ (29) 

The observed scaling of luminosity with temperature is then obtained if the 
core radius mildly decreases with increasing luminosity. A nice theoretical 
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explanation for such a scaling was first discussed by Kaiser (1991). Kaiser 
suggested that the gas in all cluster cores might have the same entropy 
as would be the case if the gas from which clusters form was preheated 
to a certain temperature. If the gas in cluster cores was then adiabatically 
compressed during the collapse of clusters, the core density should scale with 
temperature as ~ for an adiabatic index of 5/3, which translates 
into Cc ~ (Evrard & Henry 1993). This gives a scaling of luminosity 

L ~ 2’1i/ 4 (,Qj,g radius 

Tc ~ (30) 

The optical depth should then scale as 

r ~ Ue re ~ r ~ r®/® ~ (31) 
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Table 1: The estimated rms error of the x/y/z-component (galactic coordi¬ 
nates) of the bulk flow velocity of the XBACs cluster sample with respect 
to the cosmic microwave background for different cosmological models and 
two beam sizes. 





beam 

4’ 

size 

8’ 



f^bar = 0.01 

200/193/130 kms-i 

304/307/202 kms-^ 

^tot 

= 0.1 

f^bar = 0.03 

207/200/135 kms-i 

321/324/213 kms-^ 



ffbar = 0.06 

212/204/138 kms-i 

336/339/223 kms-^ 



f^bar = 0.01 

134/129/87kms-i 

211/209/138 kms-i 

^tot 

= 0.3 

f^bar = 0.03 

159/153/103 kms-i 

240/239/158 kms-^ 



ffbar = 0.06 

166/160/108 kms-i 

247/246/163 kms-^ 



f^bar = 0.01 

54/52/35 kms-i 

97/96/64 kms-i 

^tot 

= 1.0 

f^bar = 0.03 

68/65/44 kms-i 

137/134/90 kms-i 



ffbar = 0.06 

100/96/65 kms-i 

165/163/108 kms-i 
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Figure captions 


Figure 1: The la error for the determination of the peculiar velocity as 
a function of the core radius of the cluster using a “naive” hlter function 
(^(x) oc t'(x)) for three different cosmological models. The inset shows the 
angular power spectrum Ci of temperature fluctuations for the three models 
(plotted as usual as £(£ + 1)C£ against multipole £). The solid curve is for a 
universe of critical density (fltot = I), a Hubble constant of 50kms“^ Mpc~^ 
and a baryonic mass fraction of 6% (flbar = 0.06). The dotted and dashed 
curves correspond to fftot = 0.1. For the dashed curve, a non-standard reion¬ 
ization history is assumed (corresponding to an optical depth to Thomson 
scattering of order unity). The model for the gas distribution is an axisym- 
metric /3-model with core radius (T(r) oc [1 -|- (r/cc)^] [3 = 0.75). 

The results are scaled to an optical depth of the cluster of r = 0.02 and 
scale with 


Figure 2: 

a) Typical hlter function for an axisymmetric cluster in arbitrary units. The 
cluster prohle is shown for comparison. The angular scales are the same. 

b) same as a) for a non-axisymmetric cluster. The cluster is shown at the 
top, the corresponding hlter function at the bottom. 
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Figure 3: 

a) The la error in the determination of the peculiar velocity as a function 
of the core radius of the cluster using an axisymmetric “optimal” hlter func¬ 
tion (t/’(k) oc 'f'(k)/P(k)) for our hducial standard cosmological scenario 

ifg = 50kms“^ Mpc~^, and ffbar = 0.06). The beam size varies 
as indicated in the plot. The pixel noise is hxed and corresponds to 7 /uK in 
the 4’ beam . The assumed gas distribution and central optical depth are 
the same as in Fig. 1. The dashed curve shows the results for the “naive 
hlter” and a beam size of 8’ for comparison. 

b) Same as a) but for different values of the pixel noise. The pixel noise per 
beam varies as indicated in the plot. The assumed beam size is 4 arcminutes 
(FWHM). 

c) Same as a) but for a different baryonic fraction and two different beam 
sizes. Thick curve are for a beam size of 8’ and thin curves are for a beam 
size of 4’. the baryonic fraction ffbar varies as indicated in the plot. The 
inset shows the angular power spectra of temperature huctuations for the 
three cosmological models. 


Figure 4: 

a) The la error for the determination of the peculiar velocity as a function 
of the core radius of the cluster. The solid lines are for a non-axisymmetric 
optimal hlter function for our hducial standard cosmological scenario (fltot = 
1, Pq = 50kms“^ Mpc~^, and ffbar = 0.06). A scaled surface brightness 
map of Coma was used to construct the hlter function. The beam size 
(FWHM) varies as indicated in the plot. The pixel pixel noise is hxed and 
corresponds to 7 /uK in the 4’ beam. The dashed curves show the results 
for an axisymmetric hlter function for comparison (using the appropriately 
scaled best htting /3-model of Coma to construct the hlter function). 

b) Same as a) but for different values of the pixel noise. The pixel noise per 
beam varies as indicated in the plot. The assumed beam size is 4’ (FWHM). 


Figure 5: The optical depth and core radius of the XBACs cluster sample 
(Ebeling 1993) derived by using simple scaling laws, r ~ ?’c ~ 
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Figure 6: The expected la errors of peculiar velocity measurements for the 
XBACs cluster sample (Ebeling 1993) for our hducial standard cosmological 
model (fltot = 17 - 0^0 = 50kms“^ Mpc~^, ffbar = 0.06) and a beam size of 
8’ and 4’, respectively. The assumed optical depths and core radii are as in 
Fig 5. An axisymmetric optimal hlter was applied using a /3-model for the 
cluster gas distribution. The pixel noise is hxed and corresponds to 7 /uK in 
the 4’ beam. 


Figure 7: The distribution of expected la errors of peculiar velocity mea¬ 
surements for the XBACs cluster sample (Ebeling 1993) for different cos¬ 
mological models, fftot and ffbar vary as indicated in the plot. The assumed 
optical depth and core radius are as in Fig 5. An axisymmetric optimal 
hlter was applied using a /3-model for the cluster gas distribution. The pixel 
noise is hxed and corresponds to 7 /uK in the 4’ beam. 
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Figure 2b 


26 


L(T/u.u;d; - la^v )/kms 



Figure 3a 




L^T/u.u^; - iCT(,v i/Kms 



Figure 3b 












Figure 4a 


30 




200 /i,K 

X\ 

70 _ 

axisymmetric filter 

non —axisymmetric filter 

20 /LiKr^ir'=^'^-C --. ^ 


\ -' 

- 

-- 


—- 

__ 


2 /xK 


pixel 


noise 



[®core/arcmin] 


Figure 4b 








§°e o e 
*o O 

)Oo 

3 X 8’ beam 

o 4’ beam 


10 

[®core/arcmin] 


Figure 6 









